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Abstract 



^ ■ We have discovered that the eauere invariant observables of matrix mod- 

els invariant under U(A^) form a Lie algebra, in the planar large-A^ limit. 
These models include Quantum Chromodynamics and the M(atrix)-Theory 
Q^\ of strings. We study here the gauge invariant states corresponding to open 

strings ('mesons'). We find that the algebra is an extension of a remarkable 
new Lie algebra Va by a product of more well-known algebras such as gl+oo 
. and the Cuntz algebra. Va appears to be a generalization of the Lie algebra 

rS I of vector fields on the circle to non-commutative geometry. We also use a 

c^ I representation of our Lie algebra to establish an isomorphism between cer- 

tain matrix models (those that preserve 'gluon number') and open quantum 
spin chains. Using known results from quantum spin chains, we are able to 
identify some exactly solvable matrix models. Finally, the Hamiltonian of a 
dimensionally reduced QCD model is expressed explicitly as an element of 
our Lie algebra. 

PACS: n.25.Hf, ll.15.Pg, 02.20. Sv, 75.10. Jm. 

Keywords: open strings, matrix models, gauge theories, current algebra, 

quantum spin chains, large- A^ limit. 



1 Introduction 

One of the major problems in physics is to understand Yang-Mills theories 
such as Quantum ChromoDynamics (QCD). It is natural to look for a clas- 
sical approximation to this difficult problem. The obvious classical limit is 
to let Planck's constant ^ go to 0. Unfortunately, this limit is ill-suited for 
describing long-distance phenomena such as quark confinement. Another, 
less obvious, classical limit is to let the number of colors A^ go to oo. This 
large- A^ limit smooths out quantum fluctuations of gauge-invariant observ- 
ables, and is believed to manifest not only asymptotic freedom but also quark 
confinement and so should describe long-distance interactions among quarks 
better. 

We obtained the Poisson algebra of observables in this large- A^ limit in 
a previous paper |l[]. It was highly nonlinear, refiecting the dynamical fact 
that the large-A^ limit of QCD is a nonlinear classical theory. A further 
approximation is needed. This is the theory of small oscillations around the 
ground state, which can also be understood as the sum of planar diagrams 
in perturbation theory. In this paper we will study this 'planar large-A^ 
limit' of matrix models such as QCD. We will not use directly the results 
of Ref. |I[]; instead we will give a direct argument starting with the canonical 
(anti)-commutation relations. 

In a previous Letter 0, we reported on a Lie algebra describing the 
symmetry of closed string states like glueballs. The present paper serves as 
a mathematical foundation for the results in that Letter. Here we will deal 
exclusively with open string states. Building upon the mathematical results 
here, we will turn to closed string states, the mathematics of which turns out 
to be more tricky, and give a detailed exposition of the results in that Letter 
in the next paper p. 

We will not study the problems of Yang-Mills theories directly here; in- 
stead we will look at some matrix models for Yang-Mills theories. In this 
sense the present work fits into the program of studying Universal Yang-Mills 
Theories |^. The case with fermionic matter fields in Universal Yang-Mills 
theory was studied in Ref.[^. In our present language Universal Yang-Mills 
Theory corresponds to the case A = 1 with just one matrix degree of free- 
dom. This thus generalizes the considerations of Ref.]^ to the case where 
there are several matrices modeling gluons. For an earlier attempt to consider 
multi-matrix models from this point of view, see Ref.fl. 



It is well known |^, ^ that in the planar limit mesons and glueballs are 
stable, non-interacting particles. Gauge-invariant observables, such as the 
Hamiltonian, act as linear operators on these states in the light-cone for- 
malism in which (unlike in some other quantization schemes) the ground 
state is the same as that in free theory — we do not need to worry about 
non-perturbative zero mode contributions. We find that the gauge-invariant 
observables form a Lie algebra in this limit. Just as the representation the- 
ories of Kac-Moody and Virasoro algebras allow us to solve conformal field 
theories, we believe that the study of our Lie algebra will yield a solution 
to Yang-Mills theory in the planar limit. Indeed any attempt to understand 
the spectrum of QCD analytically has to begin with a solution of the planar 
limit. This means that we need to identify the symmetries peculiar to this 
limit and learn to exploit them to the fullest. 

The idea that a dynamical symmetry will help solve for the spectrum of a 
quantum theory is familiar from other problems of theoretical physics. Per- 
haps the oldest examples are the hydrogen atom and the harmonic oscillator. 
In the former, the dynamical symmetry is SO (4) and the hamiltonian can 
be written as a rational function of the generators of this Lie algebra. The 
representation theory of SO (4) then yields the spectrum. (This was Pauli's 
solution of the hydrogen atom, done even before the Schrodinger equation 
was discovered.) In the case of the harmonic socillator the symmetry is the 
symplectic Lie algebra sp(2), which is the algebra of derivations of the canon- 
ical commutation relations. The Hamiltonian is an element of sp(2). The 
familiar solution of the harmonic oscillator in terms of creation-annihilation 
operators has exactly this meaning. If we understand the dynamical symme- 
try of planar QCD equally well, we would be able to solve for its spectrum 
by a similar method. We will give an example of this in a simplified model 
of QCD, in a variational approximation. 

Conformal field theory provides more recent examples of the solution of 
a quantum theory by Lie algebraic methods |T^. The underlying symmetry 
of the Wess-Zumino-Witten model (for example) is the semi-direct product 
of the Kac-Moody algebra by the Virasoro algebra; the latter acts as deriva- 
tions on the former. The Hamiltonian is an element of the Virasoro algebra. 
The partition function of the WZW model can then be obatined from the 
Kac-Weyl character formula: this determines the eignvalues of the Hamilto- 
nian. The corresponding wavefunctions can be understood as sections of the 
determinant line on coset spaces of the Loop Group [p |. 



We are at the moment a long way from such a complete understanding 
of planar QCD. In this paper we take a crucial step in this direction: we 
identify its symmetry algebra. Somewhat surprisingly we find that it has 
a structure reminiscent of conformal field theory. Indeed it turns out that 
planar QCD is a non-commutative generalization of conformal field theory. 



The analogue of Kac-Moody algebra is the Cuntz algebra ||12| , p!3[| : it is the 
current algebra of planar QCDQ. 

There is even a natural chiral structure for the current algebra we find. 
There are two commuting (up to extensions) copies of the Cuntz algebra 
that act on the left end (the end containing an antiquark) and right end 
(the one with a quark) of a meson state. We discover the analogue of the 
Virasoro algebra [|I^] which we will call the centrix algebra. It consists of 
operators that act on the gluons in the meson state. It acts as derivations 
on the Cuntz algebra, exatly like the action of the Virasoro algebra on the 
Kac-Moody algebra. Indeed, in the special case of one degree of freedom 
for the gluons, our algebras just reduce to the Kac-Moody- Virasoro algebra. 
Even this special case is useful: it determines the spectrum of some models 
for QCD within a variational approximation. 

We discover that the symmetry algebra of planar QCD has ideals iso- 
morphic to the algebra of finite rank matrices. This shows that the essential 
spectrum of the planar QCD hamiltonian is determined by the quotient al- 
gebra. Thus the essential spectrum may be easier to determine; this is also 
suggested by the analogy with the theory of Toeplitz operators. Thus our 
mathematical study suggests some strategies for the kind of physical approx- 
imations we should make in solving planar QCD. 

It has been suspected at least since 't Hooft's paper [0 that planar Yang- 



Mills theory is closely related to a string theory. Indeed, recent works [|T7|, |18 
indicates that M-theory, the latest incarnation of string theory, is a Matrix 
model. Thus our Lie algebra is also the symmetry algebra of M-theory. A 
special case of our Lie algebra is just the Virasoro algebra, suggesting that 
our theory is indeed a generalization of string theory. 

We will find that the Lie algebra of observables has several ideals iso- 
morphic to well-known algebras such as gl+oo and Cuntz algebras and their 



^Previously, Gopakumar and Gross found that the Cuntz algebra was relevant to the 
master field [|l4| . Even in this special case, our discussion will be more complete and direct, 
being in the Hamiltonian rather than the perturbative formalism. Another previous work 
on relating the Cuntz algebra and the large- A'^ limit was done by Turgut p5[. 



products. After we quotient by these ideals we will find that what is left 
over is a remarkable new Lie algebra Va. When A = 1 this is just the Lie 
algebra of (polynomial ) vector fields on a circle, the Virasoro algebra, up 
to an extension. In the general case it is, roughly speaking, the Lie alge- 
bra of outer deruvations of the Cuntz algebra. Now, the Cuntz algebra is a 
non-commutative generalization of the algebra of functions on a circle; and 
derivations are non-commutative generalizations of vector fields. Thus we can 
think of Va as the algebra of vector fields on a non-commutative manifold, 
whose algebra of functions is Cuntz algebra. This interpretation is partic- 
ularly striking in the context of M-theory and its connection to quantum 
theories of gravity. 

If we have found an essential new symmetry algebra, it ought to help 
us solve exactly some matrix models. We are in fact able to do this by 
using our Lie algebra to prove an isomorphism between certain classes of 
matrix models (those that preserve 'gluon number') and quantum spin chain 
systems. By translating known exactly solvable spin chains into the language 
of matrix model we identify some solvable matrix models. In the case of open 
spin chains these are solutions to the Yang-Baxter and Sklyanin relations 



IS, 20, 21 



The synopsis of this paper is as follows. In Section ^, we will give a 
formulation of the physical states and the operators in the large- A^ limit. 
Then, in Section ^, we will decsribe the current algebra in the large- A^ limit. 
It has two parts — the leftix and rightix algebras of operators — that act 
on either end of a meson. In the following Section ^, we will talk about the 
algebra of the operators that act on the gluons, hence in the middle portion 
of a meson state. This will be the centrix algebra. Finally, in Section |^, we 
will piece together all the algebras we will have studied so far, together with 
the finite-dimensional matrix algebras for quarks and antiquarks, to form a 
grand Lie algebra. 

We discuss a couple of illlustrative examples of our algebra in the last 
sections. In Section |^, we will apply our algebra to study the relationship 
between quantum spin chain systems and matrix models. In Section |^ we 
will show that QCD can be formulated in terms of the grand algebra; we will 
also show to determine an approximate spectrum analytically. 

In Appendix ^, we will summarize the notations used throughout this 
article. 



2 Physical States and Operators 

We begin by introducing the basic microscopic degrees of freedom of our 
theory, which satisfy the Canonical (Anti-) Commutation Relations: 



<(^i),«i?(fc2)] = k,k25>;,i5>;,i- 


<(fci),<!(A:2)] = 0; 


«SK^i),aif(A;2)] = 0; 


[q^'{h),ql{k2)] 


^ = h.k.^l^l- 


[q^\h),q^^{k2) 


+ = 0; 


<(fci),i(^2); 


. = 0' 


'q,Aki),q^'\k2) 


^ = Kk,5>;,i- 


q^LAki),qf,2ih) 


, = 0; and 


'q^^'{k,),q^^^{k2) 


= 0. 

+ 



:l) 



(2) 



(3) 

Here, [^, -B] = AB — BA is a commutator and [A, B]+ = AB + BA an 
anti-commutator. 

The bosonic operators a, a^ create and annihilate fundamental degrees of 
freedom that we will call 'gluons'. This comes from the application of our 
theory to regularized QCD, in which case these are gluon operators. (But 



D-branes [23] or M-theory 
N label a degree of freedom 



other interpretations in terms of string bits |22 
[0 are also possible.) The indices fi,i^ = 1, 2, ■ ■ ■ 
we will call 'color' following this analogy. In the same way we will call the 
states created by g^ 'quarks' and g^ 'antiquarks'. To simplify the following 
discussion, we will mostly assume that gluons, quarks and antiquarks have 
only a finite number of distinct quantum states other than color. This can 
be done, for instance, by discretizing the momentum and other quantum 
numbers with continuous ranges of values. Let A be the possible number of 
distinct quantum states of a gluon (not counting color), and let us use the 
numbers 1, 2, . . . , A to denote these quantum states. Likewise, let us use 
A,p = 1,2, . . . ,Ap to label collectively the quantum numbers (other than 
color) of quarks or antiquarks. 

We are now ready to construct open string states, or meson states. This is 
analogous to the way closed string states, or glueball states, are constructed 
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by Thorn [Q . A meson consists of one quark, one antiquark and an arbitrary 
number of gluons. Hence, a typical colorless state of a meson can be written 
down as a linear combination of 



If' ® s^ ® (hP^ 



^-(c+l)/2-M(^^)^t..(^^^ 



a. 



i'"i+i{i'.\J 



Vi 



'{kc)ql,,{p2m, (4) 



where c = ^{K). Note that it is possible for K to be empty, which cor- 
responds to a meson state containing no gluons. Moreover, the summation 
convention for color indices is implicit in the above and following expres- 
sions. The factor of A^~('=+^)/2 is inserted so that the state has a finite norm 
in the large- A^ limit. The justification of the use of the symbol ® on the left 
hand side will be given shortly. A meson state given by Eq. (^ is depicted in 
Figs. Ka) and (b). 

Let us turn our attention to operators representing dynamical variables. 
As in models studied before (e.g., p5|, ^, |^), the only operators we will look 



at are those which can propagate single meson states to single meson states, 
in the leading order of the planar large- A^ limit. For a finite value of A^, 
they can also convert a single meson to a state with more than one meson or 
glueball. However, these terms are of order 1/A^ and so are suppresed in the 
planar large- A^ limit. This fact has been well known for a long time |2^, and 



is closely related to the planarity of Feynman diagrams of perturbation theory 



[ P^ . A non-rigorous diagrammatic proof of this fact is given in Ref . |28 . 

There are four kinds of such gauge invariant operators that survive in 
the large-A^ limit. The simplest, operators of the first kind are finite linear 
combinations of operators of the form 

^ll ® fi ® ^xt = Ar-('^+^+2)/2-tm(A^)at/;.(^,) • ■ ■al^:^^it,)qUX,) ■ 

g.,(A2X(ji)---art,(j,)g^-^(A4), (5) 

where d = #(/) and b = i^{J). In the planar limit, an operator of the first 
kind propagates a single meson state to either zero or another single meson 

state: 

Ell ® /I ® ^a! (r^ ® s"^ ® r') = K^i&' ®s'® 0^^ (6) 



It is now clear why we are using the direct product symbol ® — the 
single meson state can be regarded as a direct product of c/)^^, s^ and 0''^. 
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(a) 



(b) 



(c) 



(d) 



Figure 1: (a) A typical single meson state with 3 gluons in detail. The solid 
square at the top is the creation operator of an antiquark of Quantum State pi . 
Following the antiquark is a series of 3 gluons of Quantum States ki, k2, and 
k^ respectively. The creation operators of these gluons are represented by solid 
circles. The solid square at the bottom is the creation operators of a quark of 
Quantum State p2- Note that all the creation operators carry color indices. A 
solid line, no matter how thick it is, connecting two circles, or a circle and a 
square, is used to mean that the two corresponding creation operators share a 
color index, and this color index is being summed over. The arrow indicates 
the direction of the integer sequence K . (b) A single meson state in brief. 
The quark and aniquark are neglected. They will be consistently ignored in 
all brief diagrammatic representations. The gluon series is represented by the 
integer sequence K . (c) A typical operator of the first kind. The solid squares 
and circles are creation operators of a quark, an antiquark and gluons. The 
hollow squares are annihilation operators of an antiquark of Quantum State 
A2 and a quark of Quantum State A4, and the hollow circles are annihilation 
operators of gluons. In this particular example, there are 2 creation and 4 
annihilation operators of gluons. The creation operators are joined by thick 
lines, whereas the annihilation operators are joined by thin lines. Note that 
the sequence J is in reverse, (d) An operator of the first kind in brief. The 
thick line represents a sequence of creation operators of gluons, whereas the 
thin line represents a sequence of annihilation operators of them. J carries 
an asterisk to signify the fact that J is put in reverse. Note that the lengths 
of the two lines have no bearing on the numbers of creation or annihilation 
operators they represent. 



The set of all (p^^^s, where pi = 1,2, . . . , and Ap, form a basis of a Ap- 
dimensional vector space. The set of all 0''^'s, where again p2 = 1,2, . . . , 
and A^, form a basis of another A^t^- dimensional vector space. The operator 
of the first kind can be regarded as a direct product of the operators S^l^, 
fj and S;^^. The first operator acts as a Ap x Ap matrix on 0''^, the second 

one acts on s^ , whereas the last one acts as another Ap x Ap matrix on 0''^. 
It is therefore clear that that a meson state lies within a direct product of 
two AiT'- dimensional vector spaces (labelling the quark states) and a count- 
ably infinite-dimensional vector space spanned by all s^'s labelling the gluon 
states ( including the state containing no gluons). An operator of the first 
kind lies within the direct product gl{Ap) ® F^^gl^Ap). Here, gl{Ap) is the 
Lie algebra of the general linear group GL{Ap). Also the infinite-dimensional 
Lie algebra Fa is spanned by /j. We will see that Fa is isomorphic to the 
inductive limit gl+oo of the gl{n)^s as n — > oo. 

Operators of the second kind are finite linear combinations of operators 
of the form 



whereas an operator of the third kind can be written as a linear combination 
of operators of the form 

l^r^j^Ell ^ N-^'^'^/'ql^^ (Ai)at^^+^ ita)al^^^ (^,_i) ■ ■ ■ 

alTi^lXUlKlU2) ■ ■ ■arr^(j,)g^HA2). (8) 

Both the second and third kinds are regularized versions of the 'cur- 
rent density' operators of QCD. There are two components to the currents, 
namely the left-handed component q"'{x)q^{x) and right-handed component 
g"(x)g*(x), where a,b are spin and flavor indices. (They act on the left and 
right ends respectively of a meson.) These operators however are not well 
defined due to the divergences of quantum field theory. We must instead con- 
sider the point-split form, where the fields are evaluated at different points 
in space. But then gauge invariance forces us to insert a parallel tranport 
operator along a path P that connects the two points. Then the current 
operators become functions of an open path, i.e., string-like operators: 

rix)e'Ip''''^^iy), q'^ix)e'Ip^''^q\y). (9) 



We can now expand the field in terms of creation-anniliilation operators. 
Also, we must introduce a regularization that allows only a finite number 
of values for the momenta of quarks and gluons, in order to get a rigorous 
theory. These are precisly our operators of the first three kinds. What we 
discover is that in the large A^ limit these operators form a closed Lie algebra. 
Thus we have found the current algebra for QCD in the planar large- iV limit. 
Let us see how these operators act on a meson state, in the large N limit. 
Foe the second kind, 

S^^ ® /^ ® 1 (0^^ ® s^ ® 0^2) = 5^^ J2 6f^^'(S)s''''(g)r\ (10) 

kik2=k 

There will only be a finite number of non-zero terms in this sum (bounded 
by the number of ways of splitting K into sub-sequences), so there are no 
problems of convergence. For example, if K is shorter than J, the right hand 
side will be zero. The action of an operator of the third kind is similar except 
that it acts on the quark end: 

1 ® r J ® Ell i^'" ® s^ ® <P'") = C E ^j'^'" ® ^^'^ ® <^^' • (11) 

kik2=k 

We saw that a properly regularized version of the current operators in- 
vloves the gluon field: This means that we must also consider analoges of 
current densities that are made up entirely of gluons. In any case such terms 
will be present in the hamiltonian of QCD. They will dominate the dynamics 
of QCD in the large- A^ limit. These lead us to operators of the fourth kind, 
or gluonic operators, which are finite linear combinations of operators of the 
form: 

l®a'j®l ^ Ar-(«+^-2)/2atM.(^0at^3(,^)...^K(,^) 

<rHj6X:?a-i)---<(ji), (12) 

where a = ^{I) and 6 = #( J). Unlike the operators of the first three kinds, 
both I and J in the operators of the fourth kind must he non-empty sequences. 
In the large- A^ limit, a gluonic operator propagates a single meson state to a 
linear combination of single meson states: 

l®a^j®l{^P^®s^ ®(pp^)=^f^(d{ J2 6f^s^'^^A ®(f)f\ (13) 

\kiK2k3=k J 

10 



We see clearly from this equation that the gluonic operator acts trivially on 
the quark and antiquark. It acts on a central portion of a meson, i.e., the 
gluons lying in between the quark/antiquark pair, or on a gluon segment of 
a glueball. 

The operators defined in Eqs.(^, (H) and ([l^) are depicted in Figs. 0(c), 
(d) and |. 

As can be seen from the four kinds of operators, the essentially new 
algebraic structures are contained in the gluonic indices. In the succeeding 
sections, we are going to confine ourselves to the gluonic parts of the physical 
states and operators. The quark degrees of freedom will be put back later 
when we discuss the overall algebra, or the 'grand' algebra. 

3 Current Algebra in the Large- A^ Limit 

In this section, we are going to study the actions of the operators of the first 
three kinds on the gluon segment of a meson. This can be viewed as providing 
a representation of an associative algebra among these operators. We will see 
that the algebra of the first kind is a direct product of 2 finite-dimensional 



algebras and a well-known infinite-dimensional algebra, gl+oo p9| , |30| : it is 
thus not an essentially new object. 

Then we will turn our attention to current operators; i.e, those of the 
second and third kinds. We will see that the associative algebra among the 
current operators is a direct product of a finite-dimensional algebra and an 
infinite-dimensional associative algebra closely related to the Cuntz algebra 



1^, |T3[. Then we will focus on this infinite-dimensional algebra and study 
its structure. We will particularly be interested in the Lie algebra defined by 
the commutator in this new infinite-dimensional Lie algebra: it is a kind of 
current algebra in our theory. We will relate this current algebra with the 
Cuntz algebra |TB| and Kac-Moody algebra 



Let us consider operators of the first kind but ignoring the action on the 
quark indices; this corresponds to setting Ap = 1 temporarily. Then Eq.(|^) 
tells us that 

fjs- = Sfs^ (14) 

This equation is diagrammatically illustrated in Fig. ||(a). In particular, /g, 
where is the null sequence, is nothing but the projection operator to the 
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(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 2: (a) A typical operator of the second kind acting on the antiquark 
end in detail. There are 5 creation and 4 annihilation operators of gluons. 
The solid square is a creation operator of an antiquark, whereas the hollow 
square is an annihilation operator of it. (b) An operator of the second kind 
acting on the antiquark end in brief, (c) A typical operator of the third kind 
acting on the quark end in detail. There are 4 creation and 3 annihilation 
operators of gluons in this current operator. This time the solid square is 
a creation operator of a quark, whereas the hollow square is an annihilation 
operator of it. Compare the orientations of the arrows with those in (a). The 
choices of the orientations reflect the fact that the color indices are contracted 
differently in Eqs.l^ and (|^. (d) An operator of the third kind acting on 
the quark end in brief, (e) A typical gluonic operator in detail. There are 
2 creation and 2 annihilation operators of gluons. There are no operators 
acting on a quark or an antiquark. (f) A gluonic operator in brief. 
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gluonless state. Define the composite operator fjff by the requirement that 
its action on s^^ give fUffs^). Then it follows that 

M = ^jfL (15) 

Thus the /'s form an algebra. Moreover, it is obvious that this algebra is the 
associative algebra of finite-rank matrices on the infinite dimensional space 
spanned by the s^ . The fj are Weyl matrices in this orthonormal basis. 

Let us define a Lie algebra Fa spanned by the fj , with the usual com- 
mutator. Then 

'fjJL]=5jfl-SLff- (16) 

Hence it is clear that Fa is isomorphic to gl+oo, which is the Lie algebra, 
with the usual bracket, of all complex matrices {aij)i,j G ^+ such that the 
number of nonzero Oj/s is finite ^9). The isomorphism is given by a one-one 



correspondence between the multi-indices I and the set of natural numbers 

z+. 

Let us now turn to current operators. Ignoring the action on the (anti)- 
quark indices, we can simplify Eqs.(|l^) and ([11]) to 



/Js^ = Y. <^fs'''^and (17) 

r]s^ = E ^fs^^^, (18) 

kxk2=k 



The reader can appreciate the meanings of these equations better by looking 
at Fig. 0(b) and (c). Note that /§ = rg = 1, the identity operator. 

Let us consider the multiplication rule for the /j's. The properties of 
the algebra of the set of all rj's are completely analogous. Let us define 
the composite operator K/^ by the requirement that its action on s^'^ give 
l^jilfs^^). Then it follows that 

^M = ^Pi+ E mj.+ E ^fC^- (19) 

JiJ2=J kiK2=k 

Here we sum over all possible sequences except that J2 and K2 are required 
to be non-empty. J2 being empty and K2 being empty describe really the 
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(c) J = K2 



(d)J = K2 



Figure 3: (a) The action of an operator of the first kind on the gluon portion 
of a single meson state (Eq.^^). The dotted lines connect the line segments 
to he 'annihilated' together. The figure on the right of the arrow is the resul- 
tant meson state, (h) and (c) The actions of current operators on the gluon 
portion of a single meson state. In (b) we have the diagrammatic represen- 
tation of the R.H.S. of Eg. l\F^, whereas in (c) we have the representation 
of the R.H.S. of Eq.^^). (d) The action of a gluonic operator on a single 
meson state. 



first term on the R.H.S. and we want to make sure that this term is counted 
only once in the sums. 

Moreover, 



Thus the Vs form an algebra. This algebra is associative [^ 



this algebra has an intimate relationship with the extended Cuntz algebra 
O, |13[. Indeed, let ?i be a A-dimensional Hilbert space spanned by v^, 



v"^, . . . , v^, and J^(H) be the full Fock space ®^=o{^'"H), where ((g)°H) is 
a one-dimensional Hilbert space spanned by a unit vector Q, the 'vacuum' 
0. We will see that TiTi) is a representation space for the extended Cuntz 
algebra. There is a bijective mapping between the Fock space of physical 
states spanned by single mesons and a vector in the space J-'{TC). The state 



^For us the state fi is not exactly the vacuum, but rather is the gluonless state, which 
stih contains a quark and an antiquark. 
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s^ is mapped to the f2, and the state s^ is mapped to v^^ v'^^ ® ■ ■ ■ ® v'^". 
Define the operators a^'^ , a^\ . . . , a^^ as follows: 

a'^v^^ ® f ''^ ® ■ • • (g) ■u''" = f * ® f ''I ® f ''^ (g) . . . (g) v''" 

a'^n = v' (20) 

where i is an integer between 1 and A inclusive. The corresponding adjoint 
operators ai, 02, . . . , a\ have the following properties: 



aiv''^ ®v''^ ® ■■ -^v^" = 5'l^v^^ ® ■■ -^ v''"; and 




ttifi = 0. 


(21) 


a result, the a's and a^^s satisfy the following properties: 




Oja-'^ = 6{; and 


(22) 


y^a'^a, = 1-Pn, 


(23) 



i=l 



where Pq, is the projection operator to the vacuum Q. Thus the a's and a^'s 
are the annihilation and creation operators of the extended Cuntz algebra 
respectively. Furthermore, 

v''^ v''^ (g) ■ ■ ■ v^^ = a^'^a^^^ ■ ■ ■ a^^^n. (24) 

It is now straightforward to see that there is a one-to-one correspondence 
among the operators characterized by Eq.(^ and the operators acting on 
J-'(H). A current operator Ij corresponds to a^^^a^^t . . -a^^^aj^aj^.^ ■ ■ -ai^l^ 
corresponds to a'l'l'a^^t • • ■a*"'^; I'j corresponds to aji,aj^_^ ■ ■ ■ a-i', and l'^ to the 
identity operator. 

Define the Lie algebra between the Vs as follows: 

JlJ2 = J klK2=K 

-m- E 4jk- E ^M''- (25) 

LiL2=L iih=i 



The first three terms on the R.H.S. of Eq.(P5|) are diagrammatically repre- 
sented in Fig. ^(a) and (b). Let us call the Lie algebra defined by Eq.(p5D the 
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I -^ K ^ -71 L* 

(a) j I = k 



K 





1 

k2 

(c)j = k, 



K2 



\/ 




(b)j = kj 
/ if 

I? -■-- 
- •<2 

{ci)j,= k[ 




Figure 4: (^aj and (b) Diagrammatic representations of the defining Lie 
bracket of the leftix algebra, (a) shows the first two terms whereas (b) shows 
the third term of the R.H.S. of Eg. ([g^j are shown, (c) The Lie bracket be- 
tween l^j and f^ . Only the first term on the R.H.S. of Eq.(28) is shown. 

(d) The Lie bracket between /j and r^ . Only the first term on the R.H.S. of 
Eq.^B[) is shown. 



leftix algebra or L^- (We justify this name as an abbreviation of 'left multi- 
matrix'.) The analogous algebra spanned by rj will be called the rightix 

algebra i?A- 

Now we note the following identity: 



i=l 



lK_Ym\sM ^^M^K_ 



Thus these are nothing but ff we saw earlier: 



(26) 



if-^f-Ee 



i=i 



We can obtain from Eqs.([T^) and (|T^ the following relation: 



/ fK 



I- f 



J' J L 



E &f 



Ki rIK2 
j -I L 



KiK2=K 



LiL2=L 



I fK 
JL2- 



(27) 



(28) 



This formula is depicted in Figs. ||(c). 

As a result of Eq.(|28|), Fa forms an ideal of the algebra L\. This is 
a proper ideal as it is obvious that finite linear combinations of /j do not 
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span the whole leftix algebra. The quotient L\ = La/F^ is thus also a Lie 
algebra. Put it another way, La is the extension of La by Fa: we have the 
exact sequence of Lie algebras 

_ Fa ^ La ^ La ^ 0. (29) 

At the moment we have only this indirect construction of the Lie algebra La- 
It would be exciting to find new representations for this leftix algebra. 

How do we understand the above exact sequence in the context of the 
Cuntz algebra? The analogue of the algebra Fa is the algebra IC(H) of 
compact operators on 7i = C^; /C(7i) is just the completion of Fa in the 
topology defined by the operator norm. It follows that /C(7i) is an ideal of 
the Banach algebra generated by a, a^ . If we quotient by /C(7i), then we get 
the Cuntz algebra, with the a's and a^'s satisfying the following relations: 



ajtt 



^^ = si] and (30) 



A 



E«'^«^ = 1- (31) 



=1 



The role of the last relation is to set all /j's to zero; i.e., to quotient by Fa. 

Thus, we can regard the Lie algebra La as an extension of the Lie algebra 
associated with the Cuntz algebra by the subalgebra Fa. 

We will see that in fact the presence of Fa as a proper ideal is a generic 
feature of most of the algebras we will study. Quotienting by this ideal will 
get us the essentially new algebras we are interested in studying. But it is 
only the extension that will have interesting representations. It appears that 
the extension by gl+cx, plays a role in our theory that central extensions play 
in the theory of Kac-Moody and Virasoro algebras. Similar extensions have 



appeared in previous approaches to current algebras |^ . 

There is of course a parallel theory for Ra which is a Lie algebra spanned 
by elements of the form rj, and Ra = Ra/ Fa which is the quotient of Ra 
by Fa. Let us briefly describe it here for the sake of completeness. The Lie 
bracket expression between two r's is: 

JiJ2=J Kik2=k 

-m- E^^iJlj- E/?r- (32) 

LiL2=L Iil2=I 
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A vector in Fa can be expressed as an element in R^ as well: 



A 



/f = -f-E-:f- (33) 



1=1 



The Lie bracket between an r and an / is 

HJi]= E 5?f^''- E 4jb- (34) 

kik2=K LiL2=L 

The algebra M\ spanned by the operators Ij, rj together is also interest- 
ing. L\ and k\ are subalgebras of M\. In addition, the Lie bracket between 
an element of La and an element of R\ is 



. E . sffS- . E. ^tff£- (35) 

J1J2 = J -'l-'2 = ^ 

kik2 = k L1L2 = L 



V' ' L 



A heuristic way to see that Eq.(|35D is ture is to verify that the action of the 
R.H.S. on an arbitrary state s^^ gives lj{r^s^) — r^{l^jS^^). This, however, 
does not automatically lead to the conclusion that the Jacobi identity is 
satisfied by this definition of the Lie bracket (Eq.(|35D) because Eq.(0) and 
( |33D imply that the set of all finite linear combinations of /'s and r's is not 
linearly independent. The reader is referred to a future publication [^ for 
a rigorous proof of Eq. (|35D . A typical term in the above equation is depicted 
in Fig. 1(d). 

Therefore L\ and R\ are proper ideals of M\. In particular. Fa is a 
proper ideal of Ma, hence we have the quotient algebra M\ = M\/F\. Since 
the commutator of an / and an r is in Fa, they commute when thought of as 
elements of M\. 

These operators are regularized versions of the left- and right-handed 
current operators of QCD. We see that the chiral structure of the current 
algebra is reproduced beautifully in the multix algebra. The fact that only 
the extension of the current algebra by 5^^+00 has a representation on the 
one-meson states is also reminiscent of what happens in lower dimensional 
theories. 

Along the way we have learned an important lesson: the precisley defined 
left and right currents do not commute: their commutator is a finite rank 
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operators 


extended 
algebra 


comment 


quotient algebra 


f, 


Fa 


-^A = gi+oo 




fk _ ]k ]Kj 
■Il ~ '■l '■Li 


La 


Fa is a proper 
ideal of La. 


La = La/Fa 


fK _ ^K ^iK 
J t. ' t. 'it 


Ra 


Fa is a proper 
ideal of Ra- 


Ra = Ra/Fa 


1^1 and r^j 


Ma 


M^La®Ra 


Ma = Ma/Fa 



Table 1: Relationship among Fa, La, Ra o-nd Ma- 
tion is adopted for repeated indices in this table- 



The summation conven- 



operator instead. Our naive expectations on current algebra have to be 
revised: our extensions by Fa are just as important as the central extensions 
in the theory of Kac-Moody algebras. 

We will see that much of this structure is embedded in the algebra of 
gluonic operators on open string states. We have summarized the relationship 
among the various algebras discussed in this section in Table |l]. 

4 Lie Algebra of Gluonic Operators 

In the last section we studied the actions of the operators of the first three 
kinds, and obtained associative algebras closely related to known algebras 
such as the Cuntz algebra. It turns out that products of gluonic operators 
with each other cannot be written as finite linear combinations of these op- 
erators: they do not span an algebra under multiplication. However the 
commutator of two gluonic operators can be written as a finite linear combi- 
nation of gluonic operators. Thus gluonic operators acting on meson states 
form a Lie algebra, which we will call the centrix algebra Ea- This new Lie 
algebra we discover reduces to the Virasoro algebra (up to an extension by 
gl+oo) in the special case A = 1. Thus the centrix algebra is a sort of gen- 
eralization of the Lie algebra of vector fields on a circle to non-commutative 
geometry. 

In a previous Letter B we reported on the Lie algebra of gluonic operators 
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acting on glueball states. This is a similar but different algebra, the cyclix 
algebra. The fact that the glueball states are cyclically symmetric under 
permutations makes the cyclix Lie algebra differ in an essential way from the 
centrix algebra. 

We will see that there is a proper ideal Fa' = gl+oo of the Lie algebra S\] 
hence, there is a quotient Lie algebra E\ = X'a/Fj^. Furthermore, we will see 
that La and R\ have some elements in common with S\. This is perhaps 
a bit surprising since they originally were introduced using quark operators 
and IJ\ is the algebra of purely gluonic operators acting on meson states. 
More precisely, the generators Ij, Vj with J, J non-empty are in fact some 
linear combinations of aj. Thus we can form quotients of Sa by the algebras 
of /j's and rj's. This is the essentially new object we have discovered. 

In this section, we will only focus on the action of ctj on a segment of 
the gluonic sequence of a single meson state. We can capture the essence of 
Eq.(|T3|) by the following ^: 

a^js^= J2 S^'s^^'^\ (36) 

kiK2k3=k 

This equation can be visualized in Fig. ^d). From Eq. (|36|) , we can see that 
the set of ctj's with all possible non-empty sequences I and J's is linearly 
independent. 

Unlike the case of the algebra L\ considered in Section |^, we cannot 
define the composite operator Cja^ by the requirement that its action on 
0"*^ give (j^j{a^s^^) because in general this action cannot be written as a 



linear combination of the a's |28]. Nonetheless, the Lie bracket between two 



(j's is well-defined by the requirement that 

[-5,-f]).^ = a^(af/)-af(a5.-). (37) 



for any arbitrary sequence P. Then it can be shown [^ (by a tedious but 
straightforward calculation) that the expression of the Lie bracket is 






JiJ2=J KiK2=K 



'j'(^l 



^Of course, Uj acting on the gluonless state gives zero. 
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+ E ^"^^-fd^ E ^.Vi..+ E ^f^f' 



■J1J2 — J 

K^K2 = K 



JlJ2 = J KiK2 = J 



J J ^ J JlJ2J3=J KiK2K3=K 

K1K2 = K 
-{I ^K,J ^L). (38) 

The diagrammatic representations of the first 9 terms are given in Fig. ^ 
We will call the Lie algebra defined by Eq.(P5|) the centrix algebra T,\. 

Now let us note the following identities that follow easily from the action 



ofa5: 



-^-E^5 ^^= E ^f^^^^ 
/ Kik2=k 



\ i=l 

and 

^-5-E4)^"= E ^f^"^^- 

i=i / kiK2=k 

These are exactly the action of the opeartors Ij and rj we defined previously. 
Thus we have, 

ij = a^-E^a;and (39) 

r'j = ^^-E4- (40) 

Of course we must require that / and J be non-empty for this to be true. 

The operators Ij^s with non-empty sequences / and J form a subalgebra 
L'j^ of La, the operators Tj's form another subalgebra R'j^ of i?A, and the 
operators /j's and r j's together form yet another subalgebra M^ of Ma- What 
we have just seen is that L'^, R'j^ and Mj^ can also be viewed as subalgebras 
of X'a; in fact they are even ideals of X'a. This will be evident from the (more 
general) commutation relations between aj and Z^, etc. given below. 

Analogously, we can define a vector space F^^ spanned by fj with non- 
empty I and J. By the same counting argument as before, this is also 
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/ J^ 



/ L* 



(a) J = K 



(b) J, = K 




K2 





I 



(c)J = Ki 




\n 




KY A L* — *- I 



(e) J I = K 




'2 




(g)Jl = K, 




K 




J* K2 Y ^ L* 

(i)J = K, 







J?.-' 




^1 / \ 



K2 




K-, \ / L* 




(d)J2 = Kj 

J* Y 4> L* 

(f)J = K, 



|\"j! K 





/ J*, 





(h) J, = K 



If 

. L* 



Figure 5: Diagrammatic representations of the first 9 terms on the R.H.S. of 
Eq.^. 
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isomorphic to gl+oo- We have already seen that : 

/f = /f-E^f/' (41) 

and 

/f = rf-E^if- (42) 

4 = 1 

Thus we have 

/:^ = ^5-E4-E^^i+E^S- (43) 

j=i i=i «,i=i 

Thus Fj^ is an ideal of X'a as well. 

Now we give the commutation relations between the centrix and multix 
operators: 



-'j,if\=m+ E mj.+ E ^fr 

JiJ2=J KiK2=K 

+ E ^f^r+ E ^.i^s+ E ^^^^^f'""^ 

K1K2 = k 
-^i^f - E ^i^r- - E ^UIj 

hh=I LiL2=L 

- E siJjL2- E ^£Cj- E ^L^f.^Ls; (44) 

hh = I 

a^,rf]=5M+ E ^K.+ E ^f-r 

JlJ2=J _ffift:2=ii' 

+ E ^f-r+ E ^Sv-^+ E ^f-f^"^ 

K^K2 = k 

hh=I LiL2=L 
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LiL2=L 



I rK 



<yjJ. 



E ^l/Il^- E C^f.JL3; and(45) 

L1L2 = L L,L2L3=L 

hh = I 

(46) 






KiK2K3=K LiL2L3=L 

Eqs.(^4|), ( ^5l) and (|46|) are illustrated in Figs. P, |^ and || respectively. 








(a)J = K 



(b)Jj = K 







(c) J = K, 



(ll) J = Kj 




(e)J,= K2 




Figure 6: Diagrammatic representation of Eq.j^^). Only the first 6 terms on 
the R.H.S. of Eg. ^4^) are shown here. 

These prove, in particular, that the algebras L^, R'j^ and Fj^ are proper 



ideals of the centrix algebra. In addition, Eqs.(^J) and (|45D show that an ele- 
ment which does not belong to L'^ in the centrix algebra is an outer derivation 
of the algebra Za, i.e., the Lie algebra associated with the extended Cuntz 
algebra, and an element which does not belong to i?^ in the centrix alge- 
bra is an outer derivation of the algebra i?A. Let us summarize the various 
relationships among the Lie algebras in Table ^ 

It follows from the above discussion that there are various ways to form 
different quotient algebras from the centrix algebra. For instance, we can 
make the set of all cosets aj + F'^ into a quotient algebra E\. Furthermore, 
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(a)J = K 








It-'"" / \ / \ n 

L« I y ^\ k ^ ^- L» — »- I j^ ^^ 

(b)J2 = K 






(11) J = Kj 





(e)J,= Ki 



(0J = K2 



Figure 7: Diagrammatic representation of Eg. ^4^). Only the first 6 terms on 
the R.H.S. of Eg. (j^^ are shown here. 



operators 


extended 
algebra 


comment 


quotient 


algebra(s) 


/i 


F'k 


f'a = gl+oo 




fK _ iK ]Kj 
JL — I'L ~ i-Li 


L'a 


F'a is a proper ideal 
of L'a. 


L'a = 


L'a/F'a 


6 

fK „K „iK 

JL — ' L 'iL 


R'a 


F'a is a proper ideal 
of R'a- 


R'a = 


R'a/F'a 


Ij and rj 


M'a 


M' ^ L'a © R'a 


M'a = 


M'a/F'a 


^5 

n = ^5 - ^g 
-<y}!j + <y^ 


Sa 


F'a, L'a and R'a are 
proper ideals of Sa- 


Sa = 
Va = 


Sa/F'a 

sa/m; 

Sa/m; 



Table 2: Ea, its subalgebras and guotient algebras. The summation conven- 
tion is adopted for repeated indices in this table. 
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K2\' /\ L 



(d) J = Kj 



I Y A L* 



Figure 8: Diagrammatic representation of Eq.j^^). Only the first term on 
the R.H.S. of Eg. ^^) is shown here. 



within this quotient algebra, all cosets of the form / j + Fa span a proper ideal 
La', and all cosets of the form rK + Fa span another ideal R;^ . Thus we can 



extract the quotient Lie algebra 

Va = Ta/M; = Ta/Ma'. 



(47) 



In the simplest case A = 1, Va = Vi is just the algebra of vector fields on 
a circle. Indeed, now all integer sequences are just repetitions of the number 
1 a number of times. Write s^ as s*'^^\ Ij as Ijun, and aj as crJ^Jy (Note 



that r_^(j) 



7#W 



/^(j).) It can be easily seen that M[ is spanned by the cosets 



of the forms l^ + F[ and ll + F{, and Si is spanned by the cosets Z" + F{, 
ll + F[, 0"" + F[ and al + F[, where a and b run over all positive integers. It 
is a straightforward matter to verify the following Lie brackets, modulo F[: 



'a „1 



^1 7I 



= 






m 7l 

(2 - a - rf)(Tr'^+^ if rf < a, or 
(2 - a - d)a\_^^^ if a < d; 

(6 - rf)(T;J+rf_i; 



(c-l)/r^-^ 

f (c - l)/r''+i if 6 < c, 
\ (c-l)/L,+i ifc<6; 

-(c?-l)^6+d-i; and 
r -(rf-l)/i_,, if 

I -(rf-i)r 



-a+l 
-d+1 



if a < (i, or 
if (i < a. 



(48) 
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Let us make the following identifications: 



It - 


- z'~\ 


ll - 


^ ,-(b-l). 


^1 - 


dz 




^ -(6-1)^ 



-La-i] and 

Here 2; is a complex number with \z\ = 1. Note that the L here is not an 
integer sequence. Then Eg . (|4^) becomes 



zP, z'^' -- 


= 0; 


Lp, Lg - 

T r'? - 


= (P - (l)Lp+q] and 
— qz , 



(50) 

where p and g are integers which may be positive, negative or zero. Thus 
this is the extension of the Virasoro algebra without any central element by 
the functions on a unit circle on the complex plane. We also see that the 
algebra Vi is just the Lie algebra of vector fields on the circle, spanned by 
all Lp's for p being all integers. 

More generally, the generators of X'a act as derivations on the Cuntz alge- 
bra. When we quotient out S\ by the ideal M^, we are basically extracting 
the part which corresponds to 'inner derivations' (up to some extensions by 
standard algebras such as gl+oo-) Thus the vectrix algebra Va should be 
viewed as a generalization of the algebra of vector fields of a circle. It is the 
Lie algebra of vector fields on a sort of non-commutative generalization of 
the circle, which has the Cuntz algebra as its algebra of functions. It would 
be interesting to investigate this further. 



5 Grand Algebra Acting on Open Strings 

Now we piece together the algebras of current and gluonic operators to form 
a grand Lie algebra. In general the Hamiltonian of a matrix model in the 
large- A^ limit will be an element of this algebra. We are going to write down 
the Lie brackets of various operators, and we will see that this grand Lie 
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algebra is a sum, though not a direct sum, of four subalgebras, and there are 
two proper ideals in this grand Lie algebra. 

Let us begin by considering operators of the first kind. It is easy to see 
that all S^2 ® / j ® '^xl^^ span a Lie algebra Fa.a^^ which is isomoprphic to 
gliAp) ® Fa gl{Ap). Of course this in turn is just isomorphic to gl+oo, by 
an appropriate one-to-one correspondence with natural numbers. 

Next, let us move to operators of the second kind. We deduce from 
Eqs.(|10]) and ( p5D that the Lie bracket between two operators of this form is 

\ JlJ2 = J kiK2=K / 

-si^^®(m+ E sijL+ E ^?^rUi (51) 

V LiL2=i iih=i I 

Thus operators of the second kind form another subalgebra L\^\pOi the grand 
algebra. 

The Lie bracket between two operators of the third kind can be similarly 
derived from Eqs. ([Tl|) and ([321): 



\ JlJ2=J Kik2=K I 

-^A:i®f4ry+ E ^i/L+ E ^^rU-A^ (52) 

This equation shows clearly that all operators of the third kind form yet 
another subalgebra of the grand algebra. 

Let us consider the Lie bracket relations between operators of different 
kinds. We can derive the Lie bracket relations between an operator of the 
first kind, and either an operator of the second kind or third kind with the 
help of Eqs. (121) and (Bl) as follows: 
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S^^®/J®l,S^^®/f ®S^^ 



and 



s'^i^ii® >; 5f/r®2^:- 


-^A^H^^® >; 4./JX. 


kik2=K 


iiL2=i 


'l®r^j®Ell,Ell®fl^®Ell 


^ 




^xl^ll^ Y. '5f/f'®2^^ 


-^ii^ii^ >; ^i/£. 


kxk2=k 




LiL2=L 



2^: (53) 



2^^(54) 



These equations show that the subalgebra of the operators of the first kind 
form a proper ideal of the Lie algebra spanned by the operators of the first, 
second and third kinds. The Lie bracket between a typical operator of the 
second and that of the third kind is 



■^Ai 
"Aa 



V 



1,1 



,.K 






■^Aa 



V' ' L 



"A45 



(55) 



where 



'•p 



.K 



is given by the R.H.S. of Eq.(^). As a result, this Lie bracket 
gives a linear combination of operators of the first kind. 

Finally, consider operators of the fourth kind. The Lie bracket relations 
between an operator of the fourth kind, and an operator of any kind, are as 
follows: 







1 ® dj ® 1, 1 ® af ® 1 


= 1® 




®1; 


(56) 




l®a^®l,l®rf ®S^^ 


= 1® 


'I ^K 


®s^a; 


(57) 




l®a^j(g)l,Ell(g)lf (g)l 


= ^A^® 




(g) 1 ; and 


(58) 


"ic 


g) 


a^®l,S^^®/f ®S^^; 


= '^'Aa ^ 


g) 


'^'}J 


k 

L 


^^t 


(59) 



From the above four Lie brackets and Eqs.(P5|), (0), (^) and (^BD, we con- 
clude that all operators of the fourth kind form a subalgebra of the grand 
algebra, and that the Lie bracket between an operator of the fourth kind 
and an operator of the third, second or first kind produces an operator of 
the third, second or first kind respectively. Operators of these four kinds 
span the whole grand algebra of operators acting on mesons or open strings. 
Operators of the first, second and third kinds together span a proper ideal 
of the grand algebra. We have summarized these results in Table ^ 
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operators 


algebra 


comment 


^t ® / j ® ^t 


Fa,Af 


Fa,a^ = gliAp) ® Fa ® ^/(Af) 


^t ® ^' ® 1 


gliAp) ® La 




l®r^^®H^^ 


Ra ® gliAp) 




Ell ® ^j ® 1. and 
1 ® r^^ ® S^^ 


Ma,Ap 


Fa,Af is a proper ideal of Ma,Af ■ 


c^5 


Sa 




all of the above 


Ca^a^^ 


L'a,Ap and Ma,Ap are proper ideals of 



Table 3: T/ie grand algebra, its subalgebras and its ideals. 

6 Spin Chain Models 

Now that we have understood the algebraic properties of current and gluonic 
operators and their actions on open string states, let us give some physical 
models in which these operators and states show up. We will discuss two 
families of models in detail: mit/ti-matrix models integrable in the large- 
A^ limit associated with integrable quantum spin chain models with open 
boundary conditions in this section, and quantum chromodynamics in the 
large- A^ limit in the next section. 

In a previous Letter [0], we described the procedure for transcribing quan- 
tum spin chain models satisfying the periodic boundary condition into the 
corresponding multi-matrix models. The same procedure can be used to 
establish a one-to-one corespondence between quantum spin chain models 
satisfying open boundary conditions and multi-matrix models. The major 
modification is that boundary terms in the Hamiltonian of a spin chain corre- 
sponds to elements in La and Ra- Suppose a Hamiltonian of a matrix model 
is a linear combination H = J2i.j{hjaj + ajlj + bjrj) where hj, aj and bj are 
non-zero coefficients only if I and J have the same number of indices. Such 
linear combinations form a subalgebra of Ea] let us call it S^. Then there is 
an isomorphism between multi-matrix models whose Hamiltonians are in E° 
and quantum spin chains with interactions involving neighborhoods of spins, 
and with open boundary conditions. Furthermore, if the quantum spin chain 
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model is exactly solvable, so is the associated multi-matrix model, and vice 
versa. 

These matrix models preserve the number of 'gluons' or string bits. Such 
models can be used to study interesting string phenomena whenever the 
string tension becomes so large that the ground state favors as many string 
bits as possible ||3^. Since matrix models are appearing in many different 



contexts, the ideas in Ref.|3^ are probably intersting even outside string 
theory. 

Our work also suggests a deep relationship between integrable spin chain 
models and our Lie algebras: to each solution of the Yang-Baxter-Sklyanin 
relations, there is a maximal abelian subalgebra of the centrix algebra. This is 
just the subalgebra spanned by the conserved quantities of the corresponding 
spin chain: each conserved quantity is an element of our centrix algebra. 
Thus it emerges that the Lie algebras we have found are the underlying 
symmetries of many integrable models. We intend to develop these ideas in 
the future into a systematic theory. For now we content ourselves with a list 
of integrable multi-matrix models. 

• the six-vertex model and the corresponding integrable open quantum 
spin-^ chain model [^, ^ Q. A typical quantum state $ can be 



characterized by a sequence of c 2-dimensional column vectors, where c 
is the number of sites in the open chain. Spin-up and spin-down states 
at the j-th site are characterized by the j-th column vectors being 

and 

respectively. The Hamiltonian -f^xxz '^^ ^^^^ ^Pi^ chain model is 

1 





rj-spin 
11' 



XXZ — TT^ 

2sm7 



c-l 



+ i sin7(coth^_ri + coth^+r,")] (60) 

where 7 G (0, vr) and both ^_ and ^+ are arbitrary constants. In 
Eg . (|60|) , r|^, rj and rj are Pauli matrices at site j. They act on the j-th 
column vector. Two Pauli matrices at different sites (i.e., with different 
subscripts) are commuting. This is actually an integrable XXZ model 
with an open boundary condition. 
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We can paraphrase the model in terms of the states and operators of the 
matrix model introduced in Section ^ as follows. A single meson state 
corresponds to an open spin chain. We allow the existence of only one 
possible quantum state other than color for both quarks and antiquarks. 
Thus the notations for quarks and antiquarks in the expressions for 
single meson states can be suppressed. Each gluon corresponds to a 
spin. There are 2 possible quantum states (i.e., A = 2) other than color 
for a gluon. a^(l) (color indices are suppressed) corresponds to a spin- 
up state, whereas a^(2) corresponds to a spin-down state. Therefore a 
typical quantum state of a spin chain can be denoted by s^ , where K 
is an integer sequence of c numbers, each of which is either 1 or 2. 

The bulk term of the Hamiltonian -ffxxz^^ of the corresponding inte- 
grable matrix model in the large-A^ limit can be determined by the 
almost same procedure outlined in Ref. p|, except a minor change in 
the ranges of summations in some formulae there. The boundary terms 
on the left can be obtained from the correspondences t^ ^^ ll — I2, 
rf -|- irf <-> 2/2 and rj^ — irf <-^ 2l\. The boundary terms on the right 
can be determined analogously. Consequently, 



^^^ " ^— — |2((Ti2 +^21 +COS7(CTii -^12 -0^21 +Cr22 



22 N 
■; z I r; , „ -I- n ^1 -i- (:< )s ^y i r; , , — n ^^ — n^^ -\- n 

2 sin 7 

+i sin7[coth^_(/} - ll) + coth^+(r| - r^)]} . (61) 

We can further rewrite this formula in terms of the creation and anni- 
hilation operators a^ and a: 

^^^^^ = WV^ {2 (*(l)*(2)a-(l)<;(2) 



2A^ sin 7 

+ COS7 (at^Hl)ai-(l)a-(l)<;(l) - ajf (l)at-(2)a-(2)<^M 
- ai?(2X/(lX(l)<(2) + atf (2)at-(2)a-(2)a:;K2); 



+i sin 7 



cothe-(^^"^«i?(lX(l)gM3 - <f''aX-{2)a^^im,,) 



cothe+(gt^at^KlX(l)9'^^ - i*(2)<^(2)g'^3] J . (52) 
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integrable spin-\ XYZ model p5[. The Hamiltonian is 



c-l 



^xtz = E { [(1 + r)r;r|Vi + (1 - r)r/rjVi + Ar/T-^, 



i=i 






(63) 



where 



r = ksn^r], A = cnrjdiari, 

A± = , -D± = 



,C± 



2/i±(Aj 



sn^± sn?7± snri± 

and snw = sn(M; k), the Jacobi elhptic function of modulus < A; < 1. 
rj, C,± and X± are arbitrary constants. The corresponding Haniihonian 
in the matrix model is 



-f^XYZ = 2(cTi2 + CT21) + 2r((Tn + or. 



22 y 



+A(ai\^ - all - 4\ + ^22) + snr/ 



A_(/i^-/^) + 5_/^ + C_/? 



(64) 



integrable spin-1 XXZ model, or Fateev-Zamolodchikov model [0. The 
Hamiltonian is 



m 



Fz - E \ S, ■ S,+i - (S, • S,+i)2 + J(g - q-rS^S], 



spin 



i=i 



1 

2' 



'j "^j+i 



;(g - g-i)2(S|S|+i)2 - (g + g-^ - 2) [(^J^|^, + S]Sl^S^S]^, 



iz \2 



+l{q' - q-'){S'. - St) 



(65) 



where g is an arbitrary constant, and S^, S'^ and S^ are spin-1 matrices: 



10 ^ 

5- = 1 1 \,sy = ^ 

V 2 1 g ]^ Q ) V2 




and 
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S' 



We identify 





and 




in the spin-1 chain model by the numbers 1, 2 and 3 respectively in 
the matrix model. Then the corresponding Hamiltonian in the matrix 
model is 



--- = i-41 + crfl + 41 - all) + (--i + 41 + 41 


- 41) 


-2{all + all) - ia!l + a^ + 2all) 




—iq' + r')(-i1 + 41) + iQ + Q-')i4l + 41 + 41 


+ 41) 


+\{q - q-')\2al + 2a| - /} - r| - /f - r|) 





+ {q'-q-'){rl-rl + ll-ll). (66) 

integrahle spin-1 XYZ model. The Hamiltonian was worked out by 



Mezincescu and Nepomechie [^. We can use the same trick as above 
to work out the corresponding integrable matrix model in the large- iV 
limit. 



7 Quantum Chromodynamics in Two Dimen- 
sions 

The second family of models in which our symmetry algebras are applica- 
ble are QCD models in various dimensions, dimensionally reduced or not. 
We are going to explain briefly how to obtain the approximate Hamiltonian 
of a (2+l)-dimensional QCD model with 1 flavor of quarks and antiquarks 
by dimensional reduction from (2+1) to (1+1), which is actually a model 



of mesons studied by Antonuccio and Dalley |2^, and to write down the 
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Hamiltonian as an element of the grand algebra. We will give explicit for- 
mulae for the Hamiltonian ^. We will then make a different approximation 
from ref. ^^ (not assuming that the gluon number is conserved) which will 
give us a partial analytical understanding of the spectrum. 

In this section we will let the regulator A ^ oo, so that the momentum 
indices will take an infinite number of values. In order that the previous 
discussions of our algebra apply directly here, we will need to regularize the 
field theory such that the momentum variables can take only a finite number 
A of distinct values. But this is mostly a technicality, since we will only talk 
about field theories without divergences for which the limit A — i> oo should 
exist. 

Consider an SU(A^) gauge theory in (2+1) dimensions with one flavor of 
quarks of mass m. Let g be the strong coupling constant, a and /? be ordinary 
spacetime indices (a and P G {0, 1, 3}), A^ be a gauge potential and \1/ be a 
quark field in the fundamental representation of the gauge group U(A^). A^ 
is a traceless N x N Hermitian matrix fields whereas \l/ is a column vector 
of N Grassman fields. Let the covariant derivative -Dq,\1/ = da"^ + i^a^ and 
the Yang-Mills field be Fap = daAp — dpA^ + i[Aa, Aj^] as usual. Then the 
action is 



S 



d X 



4^ 



— TrF^flF"^ + i^7"D„^ - m^m 



(67) 



in the Weyl representation 



7 







and 7* 



T- 





(68) 



for i 

following expressions for the light-front momentum P~ and energy P"*": 



1 or 3. Following the same procedure in Ref.|]26|, we can get the 



P' 



dkk 






+1 






(69) 



roo roo 

P~ = dkhiv{k)l®al®l+ dkidk2dk'idki- 
Jo Jo 



^In fact our equation Eq.(|70|) do not agree with Eq.(14) for the self-energy in Ref.Efl]. 
But since these self-energy terms will anyway be absorbed into a redefinition of the mass 
of the scalar, they don't affect the answer. 
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.^3^4 



hivih, k2; ks, ki)5{ki + k2 - k^ - k^)! ® a'^^jf^*^ 
+hiviki; k2, k3, ki)5iki -k2-k3- k^) (l ® alf-'"'' ® 1 



-l(g)o- 



.fci 

k2k3k4 

oo 



/■OO /— ; • 



' ^S ® 2g 



+ / dkidk2dk^hn{ki\ k2, k^)5{ki - k2 - k^) {^-'^tl] 



''ki 



.r^^kl 



>l 



ki 



1 - ^\ ® /J^ 



-fci, 



rfcir 



J 



-1 ® r^^ ® s-;; + 1 ® r^^ ® s«- 






+1 (g)r 



'OO 



+ y^ dkidk2dk'idkAhn{ki,k2]k'i,k4)5{ki + k2 — kz — ki) 

(2£ j ® /g ® 1 + 1 ® rg S£ j) 

+/i//(A;i; A;2, /cs, A;4)<5(A;i -k2-k^- k^) (s^^j ® ^^^ ® 1+ 

+sfc^: ® /S^'^^ ® 1 + 1 ® r«^,3 ® s^-^: + 1 ® rl^'^ ® s^-^: 

+ / dkidk2dk-^dkjii{ki, k2] k^, k4)6{ki + k2 — k-^ — k^i] 
Jo 



(sfe: ® /, 



■0 ^ ■^fe: 



+ s 



k-i, 



'/:«-£;:) 



where 

/i/y(A;) 

hiv{ki,k2;k3,k4) 



' dp- 



hiv{ki;k2,k3,k4,] 



Ank Jo p{k — pY^ 

Svr \/kik2k^kl 

ik2 - ki){k4 - ks) {k^ + ki){ki + k2) 
{k2 + kif {k, - k2f 

Stt \/kik2k^kl 

{k2 + ki){k4-k3) (fc4 + A;i)(fc3-fc2) 



(k, + k,y 



{ks + k2y 



(70) 
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hn{k) 



hij{ki,k2;h, 



kd 




hn{ki]k2,k^,ki) 
hi{ki,k2]k^,ki) 



g'N 1 

Stt ^/k^ 

g'N 



ik2 - k^f 
2{k2 - ks' 



+ 



ki + k2 
1 



(^2 + ^3)2 k3 + k 



and 



1 



k. 



2n {ki 

(remark: hiv{k), /i/v(A;i, A;2; /cs, ^4) and /i/y(A;i; ^2, ^3, ^4) are closely related 
to the coefficients A, B and C in |2^.) The Roman numerals / and IV 
carried by some /I's refer to the fact that these are coefficients of operators of 
the ffist and fourth kinds respectively, whereas the Roman number II carried 
by other /I's signify that these are coefficients of operators of the second and 
third kinds. (Eq.([70D corresponds to Eq.(14) in Ref.|P6||, which shows terms 
corresponding to operators of the second kind only.) 

In the limit of a large number of gluons, we can regard each gluon as 
moving in the mean field created by the others. Hence it should be a sensible 
approximation to consider them all as occupying the same state with wave 
function u{k). In this case we effectively have a one-dimensional Hilbert 
space for the gluons, and our algebras simplify to the case A = 1. 

Suppose we choose an orthonormal basis where the ffist element is u{k). 
Then we can split our Hamiltonian as 



Po + Pi 



(71) 



where Pq is the contribution when all the indices correspond to u{k) and 
Pi is the rest of the terms. In fact. 



Pq = ^^<y\ + aal + pal + /3V^ 



(72) 



The last two terms in Pq describe the processes where a gluon decays into 
three gluons or three of them combine into one gluon. The coefficients are 
integrals over the wavefunction u{k): 



/^ 



dkhiv{k)u* {k)u{k); 



(73) 
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a = dkidk2dk^dkAhiviki,k2]k^.,k4)u*{k^)u*{k4)u{ki)u{k2)] (74) 

JO 

and 

P= dkidk2dk^dk4}iiv{ki]k2.,k'i,kj^)u*{ki)u{k2)u{k^)u{k4). (75) 

J oo 

We can then determine the spectrum of P^ for a given u{k) using the 
A = 1 special case of our commutation relations. It amounts to solving some 
recursion relations which we omit for the sake of brevity. 

The function u{k) can be determined by mean field theory. This is the 
same as picking the u{k) that minimizes the expectation value < P~ >, 
in each sector with a given average number of gluons. A moments thought 
will show that < P]~ >= 0. Hence we must find the spectrum of Pq^ ^^^ 
minimize the energy over all u{k). 

Clearly the sectors with even and odd numbers of gluons do not mix with 
each other. The spectrum of Pq in the even sector is 



A„ = /i + (1 + 2n)^{fi + a)2 - 4/32, n = 0,1,2,... (76) 

The 'principal quantum number' n is approximately the average number 
of gluons. The wave function is a linear combination with differing numbers 
of gluons: 

oo 
\iJn) = T.^nr\2r) (77) 

r=l 

where the coefficients $,nr depend on /i, a and f3. The lowest energy state has 
an exponentially decreasing probability to have several gluons, in agreement 
with previous studies. But as n grows, the wave function has contributions 
from states with different numbers of gluons. 

The linearly rising eigenvalues of P^ correspond to a mass spectrum that 
is also approximately linearly rising: 

M^ = A„P+ (78) 

where 

P+= / dkku*(k)u(k). (79) 

Jo 
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We can then determine u{k) by minimizing the energy for each fixed n. We 
will report more detailed information on the spectrum in later papers. 
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Appendix 

A On Multi-Indices 

Much of our work involves manipulating tensors carrying multiple indices. 
For the convenience of the reader, we give here a summary of the notations 
used in this paper for multi-indices. 

We will use lower case Latin letters such as i,j,ii,J2 to denote indices 
with values 1, 2, . . . , or A. Here A itself is a fixed positive integer, denoting 
the number of degrees of freedom of gluons. Often we will have to deal 
with a whole sequence of indices iii2^3 ■ ■ -ia, which we will denote by the 
corresponding uppercase letter I as the collective index of this sequence. 
The length of the sequence I will be denoted by #(/). 

Thus if 

/ = iii2. . .ia 

and 

J = jlJ2---jb, 

then we have #(/) = a and #(J) = b. The composition of these two se- 
quences will be denoted by 

IJ = i^i^ . . . i^jiJ2 . . . jb. 

In particular, we have 

Ij = ii?2 • • • iaj, (80) 
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when only a single index is added to the end. This operation of composing 
sequences is associative but not commutative: 

IJ ^ JI, I{JK) = {IJ)K = UK. (81) 

Often we will have to allow the null sequence among the range of values 
of a collective index. A collective index that is allowed to take the null 
sequence as its value will have a dot over it. Thus the possible values of / 
are 

J = 0, 1, 2, . . . A, 11, 12, . . . lA, 21, . . . , AA, . . . (82) 

while those of / do not include the empty set: 

/= 1,2, ■■■A,11,12,---1A,21,---,AA,--- (83) 

Of course the length of the null sequence is zero. It is the identity element 
of the composition law above, 

0/ = /0 = j. (84) 

In fact the inclusion of the empty sequence in the set of all sequences turns 
it into a semi-group under the above composition law. 

The equation I = J means that they have the same length a (say) and 

il =ji;i2 = J2;---ia =ja- 

In the same way we can define I = J either if they are both the empty 
sequence, or if they are equal in the above sense. 

The Kronecker delta function for integer sequences is defined as follows: 

ri _ j I if / = J; or 1 
^^={0 iil^J. j 

and similarly for dotted indices. The summation sign in an expression such 
as 

means that all possible distinct sequences, excluding the empty sequence, are 
summed over. (In all practical cases, it turns out that there is only a finite 

40 



number of J's such that Xj 7^ so there will be no convergence problems.) 
On the other hand, the summation sign over a dotted index 

means that all possible distinct sequences for /, including the empty sequence, 
are summed over. 

Often we will have to sum over all the ways of splitting a sequence into 
subsequences. For example, the summation sign on the L.H.S. of the equation 

J2 x^'Yk =T.^y'^^'^K (85) 

hh=i hh 

denotes the sum over all the ways in which a given index / can be split 
into two nonempty subsequences Ji and l2- If there is no way to split / as 
required, then the sum simply yields 0. For example, Eq.(|85D yields if / 
has only one integer. If the first subsequence is allowed to be empty in the 
sum, we would write instead 

E X^^Y'k- (86) 

hl2=I 

When we talk of an algebra spanned by operators such as /j, /j, r j or aj 0, 
the underlying vector space is that oi finite linear combinations. For example, 
in the linear combination X^/j cja^j, although /, J can take an infinite number 
of values, only a finite number of the complex numbers cj can be non-zero. 
Thus there are never any issues of convergence in the sums of interest to us: 
they are all finite sums. 
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